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ABSTRACT. A problem on propagation of waves in deformable 
shells with flowing liquid is very urgent in connection with wide use 
of liquid transportation systems in living organisms and technology. 
It is necessary to consider shell motion equations for influence of 
moving liquid in cavity on the dynamics of a shell by solving such 
kind problems. 

Nowadays a totality of such problems is a widely developed field of 
hydrodynamics. However, a number of peculiarities connected with 
taking into account viscous-elastic properties of the liquid and in- 
homogeneity of the shell material generates considerable mathemati- 
cal difficulties connected with integration of boundary value problems 
with variable coefficients. 

In the paper we consider wave flow of the liquid enclosed in de- 
formable tube. The used mathematical model is described by the 
equation of motion of incompressible viscous elastic liquid combined 
with equation of continuity and dynamics equation for a tube inho- 
mogeneous in length. It is accepted that the tube is cylindric, semi- 
infinite and rigidly fastened to the environment. At the infinity the 
tube is homogeneous. As a final result, the problem is reduced to the 
solution of Volterra type integral equation that is solved by sequential 
approximations method. Pulsating pressure is given at the end of the 
tube to determine the desired hydrodynamic functions. 
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Statement and mathematical ground of hydroelasticity problem is consid- 
ered that admits to describe small amplitude wave propagation in elastic tube 
nonhomogeneous in length (with modulus of elasticity E = E {x) and density 
Pm — Pm (2;)), containing liquid. The basis of liquid's model is the accounting 
of its viscous-elastic properties. Here, one-dimensional linearized equations are 
used. 

1. Given a rectilinear semi-infinite cylindrical thin-shelled tube of constant 
thickness h and radius R whose material is subjected to Hooke's law. In the 
considered case the system of hydroelasticity equations is of the form [1,2] 

du (x, t) 2 dw {x, t) _ , , 

dx dt ~ ' ^' 
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Pf^^ = -^{a{x,t)-p{x,t)}; (2) 



+ A.-^^l = 2^11 + ^.^^ j ; (3) 

ft. d'^w (x, i) 
- (a;)w(x,i) = (x) . (4) 

Rheological relation (3) describes sufficiently well the flow of liquid continua 
containing long "high-molecular" compounds. It should be noted that there 
exist two classes of variants of model (3). The media possessing instantaneous 
elasticity for which r = s + 1 belong to the first class. The models detecting 
viscous behavior at instantaneous loading belong to the second class. For them 
r = s. 

At the equations (1) - (4) u{x,t) is longitudinal speed of flow of Uquid, 
w {x, t) is radial displacement of tube's walls, p {x, t) is pressure, a (x, t) is 
"viscous-elastic" stress, p/ is liquid's density, r] is its dynamic viscosity coef- 
ficient, the quantities \j and 9j form relaxation and retardation spectra, re- 
spectively. 

Not loosing generality, we represent the functions E (x) and pm by means of 
the equalities E (x) = i?oo<7i {x) and p„L (x) = Pmoo92 {x) and accept that the 
functions gi (x) and 52 (x) are twice differentiable. We'll also assume that the 
tube is homogeneous at infinity. Hence we have: 



lim gi (x) = lim 32 (x) = '^■ (5) 

x—*oo a:— *oo 

At the same time we assume 

lim g[ (x) = 0; lim g'^ (x) = 0; Hm g'{ (x) = 0; lim g'^ (x) = 0, (6) 

X — fOO X — *oo X — >oo X — >oo 

where here and later on the primes mean differentiation with respect to x. 
Transforming the system (l)-(4) by formula 

Qix,t) ^TTR^u{x,t) 
we get the following expancc system of equation 

^^ + 2.i?^^=0, (7) 
ox ox 

pd_Q(^ = TTii^l- {a (X, t)-p (x, t)} , (8) 



.R'Y[[a{x,t)^X,'-^]^2,Y[ 



dt dx 

dQ{x,t) , ^ d'Q{x,t)- 



dx dtdx 



(9) 
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hE^ . . . , . d'^w{x,t) 

P {X, t) - -^-9l [x) W (X, t) = Pmoo%2 (x) . (10) 

Wc separate variables and look for the solution of system (1) - (4) in the 
following form: 

Q (x, t) = Qi (x) exp {itot) , w (x, t) = wi (x) exp {itot) , , , 

a {x, t) = (7i {x) exp [iut) , p {x, t) = p\ {x) exp [iLot) , ^ ' 

that allows to reduce the initial system of equations to the system of differential 
equations of the form 

Q\ {x) + 2'KRujiwi {x) = 0, (12) 
PfiujQx {x) = ttB? {x) - p\ {x)) , (13) 



^-^91 {x) - hJ^Pmoo92 (a;) I 



Pi (^) = { ~^9l {x) - hJ^pmoo92 {x) \ Wi {x) , (14) 



where in 



(1 + iLvXj) = a = qq + iai, (1 + iLv9j) = b = bo + ibi 

is accepted, and u is the given angular frequency. Further we introduce the 
denotation 

'2ri .b Cq , , Rhpmoo , ^ / hE^a \ 
PS a w 

whence 



G {x) = —uji- - —gi {x) + 92 (x) , \^co = j , (16) 



G' (x) = ^^g', (x) - ^g[ {X) . (17) 

Zpf UJ^ 



■ (17) 

Combining equations (12-15), after some transformations we get the follow- 
ing equation for the function ux (x) 

G {x) u'l {x) + G' {x) u[ (x) - ui (x) = 0. (18) 
We use the Liouville substitution 

y{x) = Q,{x)e^pi^^J^^dx^=Q,{x)./\G{xj\ (19) 
and reduce equation (18) to the form 
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y" {x) + / (.x) y (x) = 0, (20) 
where the invariant / (x) is determined by the relation 

Using relations (5) and (6) we establish the following limit equalities: 



lim G (x) = ^ - ^ + ^^^^ , lim G' (x) = 0. 



2r/ b cl ^ Bhp, 
PfLoia Lu'^ 2pf 

Hence we follow (21) and get the expression 



kp _ . ki 

1.2 1 h.2 * jt.2 1 1.2 



hm_7 (X) =6' = - i,^, (22) 



wherein 



_ 2ry aibp - apbi eg -R/tPmc 

""0 — 2~; — 2 1 — o "I n — 

PfU) a^ + af 2pf 

2ri aabi) - aibi 



Now, using denotation 



ki — 1,1 



g(x) = l-:^, (23) 



we reduce equation (20) to the form 

y" {x) + 5''y{x) = 5\{x)y{x). (24) 

We'll use a root at which ImS < later on, and on the potential q (x) we 
impose the integrability condition 

oo 

j \q{x)\dx < +00. (25) 



In order to construct the solutions, equation (24) should be completed with 
the following boundary conditions 

y (0) = 2/0, y ^ as X ^ oo. (26) 

The quantity yo depends on the functional rejime of the system, and the 
second condition (26) provides the boundedness of the desired solution. As 
a result, the solution of the hydroelasticity problem is reduced to a singular 
boundary value problem of Sturm-Lioville type (24) and (26) under the condition 
(25). 
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2. In equation (24) considering S^q (x) as an external source and applying 
the method of arbitrary constants variation we can reduce the solution of the 
problem to the equivalent integral equation [3] 

oo 

y {x, -S) = Ce-''- + 6 J sin S{^-x)q (0 y (^ , -5) d^, (27) 

X 

Here the constant C is determined as 

yo 



C 



fio,-s), 

and 

fix, -6) 

where by [4] a new function f {x,—6) is determined from the solution of the 
integral equation 



oo 

/ {x, -d)= e-''- +5jsm6{^-x)q (0 / (^ , S) d^. (28) 

X 

Equation (28) is a Volterra type equation and is solved by the sequential 
approximations method: 

oo 

f{x,-6) = J2S''fn{x,-S). (29) 
Therein we have of the following recurrent relations 



;o" (30) 

fn{x,-5)=Sj sm6{i-x)q{0fn-i{i,-S)d^ (n = l,2,...). 

X 

By inequality (25) from uniform convergences of sequential approximations 
by the Weierstrass test we can establish that the unique solution of integral 
equation (28) is determined by relation (29). We can directly establish that this 
solution solution of input equation (24) as well. 

Further, it is easy to determine the functions from systems (12)-(15) and 
write Ml, wi, pi, ai 

Q{x,t) = ^^^^yoe^-'F{x); (31) 
w{x,t) = -^^e'^'F'{x); (32) 
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a{x,t)=2T]yo-e''''F' (x) ; 



p {x, t) = yoe 



iu)t 



F'{x)., 



(33) 
(34) 



where 



F{x) = 



I [x, -S) 



^/\G{^\fo{0,-Sy 

To determine the quantity t/o at the end of the tube (a; = 0) we give the pulsating 



pressure 



p (0, t) = Po exp (iiot) . 



Using equality (34) this circumstance admits immediately to determine the 
quantity yo. It is of the form: 



yo = -« 



Po 



ai (0) ^ 92 (0) > F' (0) 



2i?w 



Thus, the solution of the problem is complete. In conclusion note that scries 
(29) in combination with relations (30) gives constructive representation of the 
desired solution and real parts (31)-(34) represent physical quantity. 
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